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' We calculate properties of A = 6 system using the accurate charge-dependent nucleon-nucleon 

^ ' (NN) potential at fourth order of chiral perturbation theory. By application of the ah initio no-core 

O . shell model (NCSM) and a variational calculation in the harmonic oscillator basis with basis size 

' up to l&hQ. we obtain the ®Li binding energy of 28.5(5) MeV and a converged excitation spectrum. 

Also, we calculate properties of ^"B using the same NN potential in a basis space of up to %hQ.. Our 
results are consistent with results obtained by standard accurate NN potentials and demonstrate a 
deficiency of Hamiltonians consisting of only two-body terms. At this order of chiral perturbation 
theory three-body terms appear. It is expected that inclusion of such terms in the Hamiltonian will 
improve agreement with experiment. 
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I. INTRODUCTION 
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' A new development in the theory of nuclear forces occurred when the concept of an effective field theory (EFT) 
was introduced and applied to the low-energy QCD Starting from a general Lagrangian consistent with the QCD 
symmetries, in particular the broken chiral symmetry, one can develop a systematic perturbative expansion valid at 
low energies for effective degrees of freedom, nucleons and pions. Pioneering work in this direction was performed by 
Ordonez, Ray and van Kolck ^, ^ who constructed a nucleon-nucleon (NN) potential in the coordinate space based 
(— I on the chiral perturbation theory at the third order. Similarly, Epelbaum et al. ^\ developed the first momentum 
■ space NN potential at the third order of the chiral perturbation theory. Recently, Entem and Machleidt developed 
.'^ [ an EFT based momentum space non-local NN potential at the fourth order of the chiral perturbation theory (next- 
I to-next-to-next-to-leading order, N'^LO) |^ that takes into account the charge dependence and fits the two-nucleon 
$H ' data with similar accuracy as the traditional phenomenological high-quality NN potentials like Argonne V18 or 
^ the CD-Bonn 2000 13 . 

It is an important question and a test for the EFT approach how well these new potentials will describe the nuclear 
structure of light nuclei when applied to more than two nucleons. The new accurate N^LO NN potential is a non- 
local momentum-space potential. There are several methods that could be applied to solve A = 3, 4 systems using 
this potential, most notably the Faddeev and the Faddeev-Yakubovsky techniques 0. However, the best established 
method for obtaining exact results for the A > 4 p-shell nuclei, the Green's Function Monte Carlo (GFMC) approach 
fej, 10, llj, is not applicable as it is best suited for work with local coordinate space potentials like AV18 or AV8' 
[l^. The N^LO potential can be used, however, in the ab initio no-core shell model approach (NCSM) j^. In this 
method, we perform calculations in a finite harmonic-oscillator (HO) basis using effective interactions appropriate 
to the basis size truncation that are systematically derived from the original inter-nucleon potential. The method is 
convergent to the exact solution with the basis size enlargement and/or with the degree of clustering of the effective 
interaction. Due to the properties of the HO basis it is no technical problem to use either local or non-local and either 
coordinate or momentum space potentials. 

In this paper, we apply the NCSM to obtain nuclear structure results using the N^LO NN potential for ■^H, ^He, 
6Li, 6He and i°B. With the improved shell model code ANTOINE [11 [13 we are able to reach basis sizes up to the 
16?ir2 for ^Li and up to the Shft for ^"B. At the same time, we perform variational calculations in the HO basis for ^Li 
using the bare N'^LO potential. These variational calculations allow us to make extrapolations and check our NCSM 
effective interaction results. In addition, we compare the N^LO results to those obtained by the CD-Bonn 2000 and 
the AV8'. 

It should be pointed out that already at the third order of the chiral perturbation theory, N^LO, three-body terms 
appear These terms should be consistently included in calculations for A > 2 nuclei. It is our goal to perform 
the NCSM calculations with the three-body terms included using the approach outlined in Ref. ^^I- In f^ict, work in 
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this direction is under way 18]. However, in general when only the two-nucleon interactions are employed the NCSM 
approach can achieve a higher accuracy as larger basis spaces can be reached. Therefore, it is valuable to perform 
accurate calculations using just the two-nucleon forces to assess the effects of the omitted three-nucleon terms. The 
eventual less accurate results obtained with both the two- and three-nucleon forces can then be better understood 
and extrapolated. 

We briefly review the basic features of the NCSM and describe our variational calculations in the HO basis in 
Section im In Section Hill we present our N'^LO results for ■^H and ^He. In Section Hvl we discuss our calculations 
for A = 6 using the N'^LO NN potential and compare them to results obtained by the AV8' and the CD-Bonn 2000. 
In addition, we check convergence properties of our approach using the semi-realistic Minnesota NN potential. In 
Section^ we show the NCSM predictions for ^"B with the N'^LO NN potential. Our conclusions are summarized in 
Section EH 



II. AB INITIO NO-CORE SHELL MODEL AND VARIATIONAL CALCULATIONS IN THE HO BASIS 



A detailed description of the NCSM approach was presented, e.g. in Refs. |13ll22ll23| . Here we only briefly review 
the basic features of the NCSM. Also, we describe our variational calculations in the HO basis that we apply in the 
following Sections. 

The starting Hamiltonian for our investigations is 
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where m is the nucleon mass, VNN.ij, the NN interaction with both strong and electromagnetic components. In 
the NCSM, we employ a large but finite harmonic-oscillator (HO) basis. Due to properties of the realistic nuclear 
interaction in Eq. we must derive an effective interaction appropriate for the basis truncation in order to 

reach convergence or to achieve a reasonable approximation to the exact solution. To facilitate the derivation of 
the effective interaction, we modify the Hamiltonian ^ by adding to it the center-of-mass (CM) HO Hamiltonian 



TcM + UcM, where Uq 
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This modification was first introduced by Lipkin 



[l9l |. The effect of the HO CM Hamiltonian will later be subtracted out in the final many-body calculation. Due 
to the translational invariance of the Hamiltonian the HO CM Hamiltonian has in fact no effect on the intrinsic 
properties of the system in the infinite basis space or in a finite basis space truncated by A^max as described below if 
the interaction in Eq. ^ is not altered. 

The modified Hamiltonian can be cast into the form 



Ha + HcM ^Y.'^^+Jl^^ 
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Next we divide the A-nucleon infinite HO basis space into the finite active space [P) comprising of all states of up to 
A'max HO excitations above the unperturbed ground state and the excluded space {Q = 1 — P). The basic idea of the 
NCSM approach is to apply a unitary transformation on the Hamiltonian (|2Jl, e~^H^e^ such that Qe~^ H^e^ P = 
[23. 125I I . If such a transformation is found, the effective Hamiltonian that exactly reproduces a subset of eigenstates 
of the full space Hamiltonian is given by Hcs = Pe~^ H^e^ P. This effective Hamiltonian contains up to A-body 
terms and to construct it is essentially as difficult as to solve the full problem. Therefore, we apply this basic idea on 
a sub-cluster level. In this paper we use the simplest approximation, a two-body effective interaction approximation. 
This approximation is obtained by constructing the unitary transformation for the Hamiltonian ^ applied to two 
nucleons only, i.e., 7^2 = hi + h2 + V^'^ with A in V^'^ fixed at the value corresponding to the final A-nucleon 
calculation. The two-body effective interaction is then derived with the help of the exact eigensolutions of Ti,2 as 
[^i^'^lcfi — P2 [e~'^^7Y2e2 — hi — /12] P2- Here, the definition of the two-nucleon model space projector P2 follows 
directly from the definition of the A-nucleon projector P. The two-nucleon transformation S2 is determined by the 
condition Q2e~^''H2e^'^ P2 = 0. 

The A-nucleon calculation with A > A is then performed using the multi-shell (no-core) version of the code AN- 



3 



TOINE, with the Hamiltonian in the form: 

+ /3(i^CM-^?^^^)|p, (3) 

where we subtracted the CM Hamiltonian Hcm and added the Lawson projection term f3{HcM — ^Ti^) to shift the 
spurious CM excitations. As our effective interaction is translationally invariant our physical eigenenergies do not 
depend on the choice of (3. Note that in Section ITTll we present A = 3,4 calculations performed in Jacobi coordinate 
antisymmetrized HO basis. In those calculations the CM degrees of freedom are explicitly omitted and therefore no 
fi{HcM — ^h^) term is used. 

We note that our effective interaction depends on the nuclcon number A, the HO frequency fJ, and the P-space 
basis size defined by A^max- It is translationally invariant and with iVmax oo the NCSM effective Hamiltonian 
approaches the starting bare Hamiltonian Consequently, with the basis size increase the NCSM results become 
less and less dependent on the HO frequency and converge to the exact solution. Alternatively, by increasing the 
clustering of the effective interaction for a fixed P-space size the NCSM effective Hamiltonian approaches the exact 
A-nucleon effective Hamiltonian that reproduces exactly a subset of eigenstates of the starting Hamiltonian . 

It turns out that the N^LO NN potential is softer than the phenomenological high-quality NN potentials like the 
AV18 or the CD-Bonn. As we are able to reach basis spaces of up to A^max = 18 and A^max = 16 for A = A and 
A = 6, respectively, it make sense to perform variational calculations in the HO basis for ''He and ^Li using the bare 
N'^LO NN potential with the HO frequency or the oscillator length as the only variational parameter. In this case, 
the Hamiltonian that we use is PHaP, with Ha given in Eq. iP) (plus the pP{HcM - ^hn)P term when Slater 
determinant basis is utilized, i.e., for A > A). We perform variational calculations for different basis sizes defined by 
A^max- We then extrapolate to A^^ax ~^ oo and compare to our effective interaction results. 
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III. AND ''HE RESULTS 



As a test of con verg ence of our method and also as a test of correct interfacing of our codes with the N'^LO 
NN potential code |2y| we performed calculations for ^H and ^He. We used the translationally invariant HO basis 
antisymmetrized as described in Ref. and employed the many-body effective interaction code MANYEFF [23l |. 
Our results are summarized in Tabled Our ^H binding energy, 7.85(1) MeV, that was obtained using basis spaces 
up to A'lnax — 40 agrees well with the result obtained by the Faddeev method [2(tI |. 

Our "^He binding energy calculations were performed using both bare and the two-body effective interactions in 
basis spaces up to A^^ax = 18. Our obtained binding energy, 25.36(4) MeV, is very close to that obtained by the 
Faddeev- Yakubovsky method In Fig. ^ we present the HO frequency dependence of our ''He binding energy for 
different basis sizes. Calculations with two-body effective interactions are compared to the variational calculations 
using the bare N^LO NN potential. The NCSM calculations with the two-body effective interaction are not variational 
as some terms of the transformed Hamiltonian (i.e., the higher than two-body terms) are omitted. Therefore, in the 
NCSM approach the convergence could be, depending on the HO frequency, either from above or below or even 
oscillatory with the basis size enlargement. It should be noted how the frequency dependence of our '''He binding 
energy decreases with the basis size enlargement. Also, the two-body effective interaction always improves on the bare 
interaction results in particular for the smaller spaces and lower HO frequencies. For example, the A'max — 14 {lAMl) 
binding energy obtained using the effective interaction varies within 25.11 MeV and 25.5 MeV in the range of HO 
frequencies shown in Fig. ^ UVt = 24 — 40 MeV. The bare interaction binding energy in the same basis space and the 
same frequency range varies within 23.09 MeV and 25.12 MeV. For the highest frequency that we used here. Ml = 40 
MeV, we observe that the improvements due to the two-body effective interaction become very small. We note that 
the the minimum of the variational calculation in the A^max ~ 18 basis space is at about hil. = 36 MeV. The minima 
for the smaller space calculations are at still higher frequencies, e.g., the minimum of the of the A'max = 12 calculation 
is at hfl > 40 MeV outside of the range shown in Fig. ^ At A^max — 18 our effective interaction calculations agree 
with the variational bare interaction result. 

It is instructive to investigate in addition to the binding energy also the convergence of the point-nuclcon ^He 
radius. A similar plot as in Fig. ^ is presented for the point-nucleon root mean square radius in Fig. |21 Even 
though no operator renormalization was used for the radius operator very similar conclusions can be drawn as for the 
binding energy. The calculation with the effective interaction (i.e., with wave functions obtained using the effective 
interactions) improves convergence and reduces the frequency dependence compared to a calculation with the bare 
interaction wave functions using the same basis size. With the basis enlargement the frequency dependence decreases. 
At the same time, it should be noticed that the optimal frequency for the radius convergence is not necessarily the 
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same as that for the binding energy. It is important to investigate the frequency dependence of energies and other 
observables to determine the extrapolated values and errors. We extrapolate the ''He radius with the N'^LO potential 
to be 1.515(10) fm. In Fig. Owe show the complete basis size dependence of the point-nucleon radius calculation 
using the effective interaction for several HO frequencies. The convergence with iVmax and the decrease of frequency 
dependence can be clearly seen. 

The and ^He N^LO binding energies (7.85 MeV, 25.36 MeV) are in between the Argonne V8' (7.77 MeV, 25.2 
MeV) and the CD-Bonn binding energies (8.00 MeV, 26.3 MeV) but considerably closer to the AV8' results. Note 
that the Coulomb interaction is included for all three potentials. 

IV. APPLICATION TO '^LI AND ®HE 

As mentioned earlier, our A = 6 calculations reported here were obtained using the improved multi-shell (no-core) 
version of the code ANTOINE capable to reach basis spaces up to A'max = 16 (16?ifi) for ^Li with the M-scheme basis 
dimension equal to 8 x 10^. A specific challenge of the no-core calculations compared to the traditional Ohn, valence 
nucleon calculations is the enormous number of single particle states (e.g., there are 171 nlj levels and 2280 nljm 
states for each protons and neutrons in the 16hil ^Li calculations). This implies a huge number of operators that 
need to be stored in memory. The other challenge is the huge dimension of the matrix to be diagonalized. Both these 
challenges have been addressed in the improved ANTOINE code. In particular, an efficient way of operator storage 
was found and a partitioning of Lanczos vectors was introduced that eliminates the need to store the full Lanczos 
vectors in memory. 

Before presenting our N'^LO results, let us first discuss test calculations for ^Li using the semi-realistic Minnesota 
(MN) NN potential 26] frequently used for few-body calculation benchmarks 113, H^- Our NCSM results, in particular 
the HO frequency dependence of the ground-state energy for different basis spaces ranging from Nyna^ = to A^max = 
14, are presented in Fig. ^ Clearly, for the Minnesota NN potential we observe the best case scenario of the NCSM 
convergence. Starting from A^max — 2 the convergence is uniform in the whole frequency range. The HO frequency 
dependence is getting weaker with the basis size enlargement. The differences between the successive curves decrease 
as the basis size increases. At A^max = 14 the NCSM energy at the minimum is within 30 keV of the Stochastic 
Variational Method (SVM) ^ result 36.51 MeV An Af,„ax = 16 calculation improves the NCSM and SVM 
agreement to 10 keV. Similarly, the Effective Interaction in the Hyperspherical Harmonics basis (EIHH) method 
[28Ll30| obtains a very close binding energy of 36.64(7) MeV jsj. Let us mention that a variational calculation in the 
HO basis using the bare Minnesota NN potential still misses about 1.5 MeV from the converged value in the A^max — 14 
basis space. The improvement due to the two-body effective interaction is quite dramatic even for such a large basis 
space. The Minnesota NN potential was employed as a test potential for the NCSM ^Li investigation already in Ref. 
[3^ . There, however, basis spaces up to only A^max = 10 were utilized. Also, the Minnesota NN potential results 
presented in Ref. [l^ include the Coulomb interaction which was omitted in the present calculations. 

We now turn to our ^Li N'^LO results. In Fig. [SJ we present the HO frequency dependence of the ground-state 
energy obtained by the variational calculations (dashed lines) in A'max = 8 to A^max = 16 basis spaces and by the two- 
body effective interaction NCSM calculations in basis spaces ranging from A^max = 2 to A^max = 14. The variational 
calculations exhibit large changes with A'max and clearly they are not converged at A^max = 16. The minimum in the 
largest spaces that we used develops at about Tifl — 33 MeV and the value at minimum for the A^max = 16 space 
is -24.7 MeV. The NCSM two-body effective interaction results (full lines in Fig. O show stronger HO frequency 
dependence than those obtained using the Minnesota NN potential and in fact also than those obtained using the 
CD-Bonn NN potential [s^- We can see that the minimum appears around Tifi = 12 MeV and shifts with increasing 
A^max to lower frequency. Due to this shift the successive curves intersect and we do not have a nice uniform picture as 
in the case of the Minnesota interaction. At the same time, the value at minimum does not change much starting at 
A^max = 6 as can be seen in Fig. Elwhere we plot the ground-state energy obtained at the minimum for each A^max- It 
is interesting to note that for the largest spaces, e.g. A^max = 14, a second minimum develops at around the frequency 
where the bare interaction results are at the minimum. Similarly as in the case of ^He, we can see that starting at 
certain frequency, the two-body effective interaction does not change the bare interaction results any more and the 
effective interaction curves start to follow the bare interaction curves. 

In Fig. 1^1 we also present the bare interaction variational calculation results at fixed HO frequencies at or close 
to the minimum of the A^max = 14 and A^max — 16 spaces. It turns out that these ground-state energy A'max 
dependencies can be fitted by an exponential function of the t ype ii^(A',„ax) — Eoo + ae^''^'""^ (see an analogous fit 
done for the hyperspherical harmonics basis calculation in Ref. \2&^). When we fit this formula to the last four points, 
A^max = 8 — 16, we obtain Eoo ~ —28.2 MeV. This is rather close to the energy at the minima of our two-body 
effective interaction calculations that are at about 28.5 MeV. Based on the variational bare interaction calculations 
and extrapolation and on the two-body effective interaction NCSM calculations in different spaces and using different 
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frequencies we arrive at our final ^Li N'^LO binding energy result of 28.5(5) MeV. The experimental ^Li binding energy 
is 31.995 MeV. Our calculations suggest that the N'^LO NN potential under-binds ^Li by about 3.5 MeV. 

It is interesting to compare the N'^LO results to those obtained by other NN potentials. Concerning the binding 
energy, in Fig. Hwe repeat the basis size dependence from Fig. Elof the NCSM ground-state energy at the frequency 
minima for the N^LO and show the same for the AV8', CD-Bonn 2000 and the Minnesota NN potentials. In addition, 
the AV8' GFMC result 10] and the Minnesota SVM 29J and EIHII ^31,] results are shown for comparison. As discussed 
earlier, our Minnesota results show a nice convergence and agree quite well in particular with the SVM result. Our 
AV8' binding energy is close to that obtained using the N'^LO NN potential and within about 300 keV of the GFMC 
result. Our CD-Bonn 2000 binding energy is, on the other hand, larger by almost 1 MeV compared to the N'^LO and 
the AV8'. Our present result is almost the same as that presented in Ref. [s^l using the older version of the CD-Bonn 
NN potential and the NCSM calculations in basis spaces up to only 10?ir2. Our binding energy results are tabulated 
in Tables Hll and ITTll with our conservative error estimates based on the frequency and basis size sensitivity. 

We investigated the excitation energies of the five lowest excited states of ^Li using the N^LO NN potential. As 
the binding energy minimum shifts with the change of the basis size and also as we expect that the energy of different 
excited states might have different dependence on the HO frequency than the ground state, we calculated the excitation 
energies for several HO frequencies. Due to the complexity of the calculations, the lowest four state were obtained in 
basis spaces up to A^max = 14 while we stopped at A^max = 12 for the 2+1 and the 1^0 state. The NCSM excitation 
energy dependence on the basis size is presented in Figs. |H1 IHl and ^| for the HO frequencies of hil — 8, 10 and 13 
MeV, respectively. In Fig. ^] we then show the excitation energies for the hft = 8, 10, 12 and 13 MeV HO frequencies 
in the largest basis space used. It is apparent that the convergence rate with A^max is different for different states. 
In particular, the 3+0 state and the 0+1 state converge faster in the higher frequency calculations, while the higher 
lying states converge faster in the lower frequency calculations. As can be seen from Fig. in the largest space 
used the dependence on the HO frequency is very weak. We summarize our excitation energies with the estimated 
errors in Table ^ We have, in particular, confidence in our 3+0 excitation energy that is reflected in a very small 
associated uncertainty. In the A^max = 14 space, the 3+0 excitation energy is almost frequency independent in the 
range of hft = 10 — 13 MeV. At hfl = 8 MeV, this state has a slightly higher excitation energy. However, it is seen 
from Fig. |Slthat its energy is still decreasing with the basis size enlargement. On the other hand, this state is fairly 
stable and converged with A'niax in the hfl = 10 — 13 MeV calculations. 

In general, our calculated levels are in the correct order when compared to experiment for the lowest four states. 
Concerning the 2+1 and the l2'0 state, those are obtained at almost the same energy and the uncertainty of our 
results is too large to make a definitive prediction as to their ordering. When comparing our results to experiment 
we notice a striking discrepancy for the 3+0 state that calculated excitation energy is about 0.7 MeV higher than the 
experimental one. Taking into account our confidence in convergence of this state, this is a serious discrepancy in 
the prediction of the N^LO NN potential. The remaining states, when considering our estimated errors, are in fairly 
good agreement with experiment. We note, however, a significant underestimation of the splitting between the 3+0 
and the 2+0 states. This splitting is influenced by the spin-orbit interaction strength. 

In Table lllll we compare our N'^LO excitation energy results with our NCSM results obtained by the CD-Bonn 
2000 and the AV8' NN potentials. In addition, we also present the GFMC ""Li results obtained using the AV8'. We 
note that we have not performed as detailed excitation energy calculations using the CD-Bonn 2000 and the AV8' 
as we did for the N'^LO. The presented NCSM CD-Bonn 2000 and the AV8' excitation energies were obtained in the 
A^max = 14 basis space using the HO frequency of hft = 12 MeV and Ti^ = 11 MeV, respectively. We expect, however, 
that the error estimates that we found for our N'^LO excitation energy results are also applicable to our tabulated CD- 
Bonn 2000 and AV8' excitation energy calculations. It is apparent that when considering the estimated uncertainties 
our NCSM calculations cannot discriminate between the three realistic NN potentials as to their predictions of the 
excitation energies. The CD-Bonn gives a slightly lower excitation energy of the 3+0 state and a slightly larger 
splitting between the 3+0 and the 2+0 states. However, the differences are smaller than the difference between our 
NCSM AV8' and the GFMC AV8' results. In general, the NCSM and the GFMC results using the AV8' are in a 
reasonable agreement. Although we have an about 200 keV difference in the 3+0 excitation energy, the difference 
with respect to experiment is about 800 keV for the NCSM and about 1 MeV for the GFMC. This further supports 
our statement about the failure of the accurate NN potentials to predict the position of this lowest excited state of 
6Li. 

The ^Li with the CD-Bonn and the AV8' was investigated within the NCSM in Refs. [s^ls^ using the basis spaces 
up to A^max = 10 with the two-body effective interaction and A'max = 6 with the three-body effective interaction. 
The results that we present here are consistent with the previously published results. Only the 2+0 excitation energy 
reported in this work is better converged and smaller. This is due to the larger basis space, i.e. A^^ax — 14, employed 
here. We note that there is approximately a 500 keV difference in the excitation energy of the 0+1 state obtained 
in the NCSM and the GFMC using the AV8'. It is interesting to point out that the NCSM calculations using the 
three-body effective interaction in the A^max — 6 space predict this 0+ 1 state at a higher excitation energy closer to the 
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GFMC result than the NCSM calculations that employ just the two-body effective interaction in similar basis spaces 
[ssf . This is contrary to the T = states excitation energy results that are very similar using both the two-body or the 
three-body effective interaction [s^. The difference for the 0+1 state in the two NCSM approximations needs further 
investigation. Ultimately, with the basis size enlargement the results in both approximations must by construction 
converge to the same value. 

In Table im we also present electromagnetic properties of ^Li, binding energy of ^He and the Gamow- Teller ^He^^Li 
B(GT) value obtained in our NCSM calculations using the N^LO NN potential. The ^He binding energy is consistent 
with our ^Li results and under-binds experiment by about 3 MeV. The NCSM calculated electromagnetic properties, 
obtained using bare charges, are in a reasonable agreement with experiment. We note that while the B(M1) result is 
fairly stable with the basis size change the B(E2) values increase in general as the basis size is increased. Consequently, 
our B(M1) values reported here are almost the same as those published in Ref. [s^l but our present B(E2) values are 
substantially larger. The present results were obtained using N^ji^x = 14 while those of Ref. using TVmax = 10. 
We note that our B(GT) value for the ^He— >^Li ground-state to ground-state transition over-predicts experiment. 
Recently, Schiavilla and Wiringa found that the AV18/Urbana-IX interaction also over-predicts the ^He^^Li B(GT) 
value III iB,ct, their obtained Gamow- Teller matrix elements (= ■\/B(GT)), 2.254(5) and 2.246(10) using two 

types of wave functions, compare well with our result 2.28 obtained using the N'^LO NN potential. 

In order to judge the convergence of other observables than just energies, in Fig. 1121 we present the HO frequency 
dependence of the ^Li ground-state quadrupole moment for different basis sizes. Similar conclusions as for the '*He 
radius calculations can be drawn here. With the basis size increase the frequency dependence decreases. At the same 
time, the fastest convergence is not necessarily obtained for the same frequency as for the binding energy. A speed up 
of convergence for the electromagnetic quadrupole observables should be achieved by using effective operators. This 
has not been done here, but work on this problem is currently under way l35|. We conclude that the quadrupole 
moment of ^Li with the N'^LO interaction is -0.08(2) e fm^. 

To summarize this section, apart from more binding obtained using the CD-Bonn 2000 NN potential, we can see 
little differences in predictions of the A = 6 nuclei properties by the three accurate NN potentials, the N^LO, the 
CD-Bonn 2000 and the AV8'. Let us just note that we also investigated the 2?ir2-dominated (intruder) states in 
^He (not observed in experiment) using the three NN potentials. Such states observed in heavier p-shell nuclei are 
in general predicted by the NCSM at too high energies and converge much more slowly than the p-shell dominated 
states 0, . Using the same basis size and the same HO frequency, these states appear at lower excitation energy 
when the N^LO NN potential is used compared to the CD-Bonn 2000 and the AV8'. This does not necessarily imply 
that the N'^LO NN potential will predict the 2?ii7-dominated states at lower energies but perhaps that the NCSM 
convergence of such states is faster when the N^LO NN potential is used. 



V. APPLICATION TO 



In the previous section we pointed out that the N'^LO NN potential fails to predict correct excitation energy of 
the 3+0 state in ^Li by about 700 keV. It has been realized that this problem which is common to other accurate 
NN potentials is magnified in ^'^B. The ground state of ^°B is 3+0 and the 1+0 state is the ^°B first excited state 
at 0.72 MeV. Both the NCSM calculations using the CD-Bonn and the Argonne NN potentials and the GFMC 
^"B calculations using the Argonne NN potentials show that these accurate NN potentials predict incorrectly the 1+0 
state as the ground state and the 3+0 state as an excited state at higher than 1 MeV of excitation energy pTLl3a . l3a |. 
This is a clear indication for the need of a three-nucleon force not only to fix the under-binding problem but also to 
correct, at least in some nuclei like ^'^B, the level ordering. In fact, it has already been shown that, e.g., for the AV8' 
NN potential the ^°B ground-state spin problem is resolved by including certain types of the three-nucleon forces like 
Illinois [Till or the Tucson-Melbourne TM'(99) [13, but it is not resolved by including the Urbana IX ^| . This type 
of sensitivity suggests that nuclear structure calculations for light nuclei could be used to help determine the form and 
the parametrization of realistic three-nucleon forces that are unlike the two-nucleon interactions not well established. 

It is quite interesting to perform the ^°B calculations using the N^LO NN potential. Our NCSM results performed 
in basis spaces up to iVniax = 8 are presented in Figs. 1131 and [TH and Table HVl Unlike in the case of ^Li, we were 
not able to reach as large basis spaces that are needed for the excitation energy or the binding energy convergence. 
Still, our calculations are sufficient to resolve the ground-state spin issue in i°B with the N^LO NN potential. In Fig. 
1131 we present the dependence of the ^"^B 1+0 and 3+0 state energy on the HO frequency for different model spaces. 
In the NCSM, one typically determines the optimal HO frequency for the ground state and uses the same frequency 
to describe the excited states. In general, however, one should also investigate the the frequency dependence for the 
excited states in order to determine their energy p33 |. This is in particular important in the ^"^B investigation when 
we actually want to find out which state is the ground state for a given NN potential. In the inset of Fig. ^1 the 1+0 
and 3+0 state energies at their respective HO frequency minima are plotted as a function of A^max- We can see that 
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in the whole frequency range that we show the 1+0 state is always the ground state. For the highest frequency shown 
in the Figure, hH. — 16 MeV, the excitation energy of the 3"''0 state becomes quite small. However, we observe that 
with the basis size increase the 3~''0 excitation energy is actually increasing. In Fig. 1141 we then show the excitation 
energies of the lowest 5 states calculated in the Ohfl — 8?ifi basis spaces using a fixed HO frequency of hfl = 12 
MeV. At this frequency, the ground-state energy is at the minimum in the Shfl space. Clearly, the excitation spectrum 
is not converged as it was the case in our ^Li calculations. However, there is definitely a trend to convergence. The 
differences between the successive calculations become smaller and smaller with increasing A^max- We can conclude 
with confidence that the 3"'"0 state will remain the excited state with further increase of Amax and most likely will 
appear at an excitation energy of more than 1 MeV. Based on our NCSM calculations, we conclude that the N'^LO 
NN potential fails to produce a correct ground-state spin for ^"B. 

We mentioned earlier that already at the N^LO order of the chiral perturbation theory the three-body terms 
appear. Those terms were not used in this work. However, work on their inclusion in the NCSM is under way jl8j |. It 
is reasonable to expect that these three-body terms will resolve the ^°B ground state problem. Structure of some of 
the EFT three-body terms in fact coincides with, e.g., the Tucson- Melbourne TM'(99) three-nucleon interaction [s^ 
that was already shown to resolve the ^"B problem when combined with the AV8' NN potential. 

VI. CONCLUSIONS 

We calculated the properties of ^^H, "'He, ^Li, ^He and ^"B using the accurate charge-dependent NN potential at 
fourth order of chiral perturbation theory. We applied the ab initio no-core shell model and a variational calculation 
in the harmonic oscillator basis. For ^Li, we were able to reach the basis size of up to IGhQ, using an improved 
multi-shell (no-core) version of the shell model code ANTOINE. We obtained the ^Li binding energy of 28.5(5) MeV 
and a converged excitation spectrum. Properties of ^°B were obtained using the same NN potential in basis spaces of 
up to 8?ir2. Our results are consistent with results obtained by standard accurate NN potentials and demonstrate a 
deficiency of Hamiltonians consisting of only two-body terms. In addition to under binding compared to experiment 
for all investigated nuclei, we found the ^Li 3"'"0 excitation energy over-predicted by about 700 keV compared to 
experiment. For ^"^B we found an incorrect ground state, 1~''0, contrary to experimental 3"'"0. 

We anticipate that most of these problems will be resolved by including three-body terms that appear already at 
the third order of the chiral perturbation theory. A work on including these terms in the NCSM is under way |ia |. 

Finally, let us note that recently a new accurate non-local NN potential has been constructed that fits not only the 
two-nucleon data but also the binding energies of "^H and ■^He without any need to introduce a three-body interaction 
|38j |. As our calculations show, the deficiencies of two-nucleon interactions are not limited only to the under binding 
but also to the nuclear structure issues most likely linked to an insufficient spin-orbit interaction strength. It will be 
very interesting to perform nuclear structure calculations using the new non-local NN potential of Ref. to test 
its predictions for the binding energies and also for the spin-orbit interaction strength sensitive observables. 
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FIG. 1: (Color online) Harmonic-oscillator frequency dependence of the *He ground-state energy obtained in 12hQ-lShQ 
(A'^max = 12 — 18) basis spaces using the bare N^LO NN potential (dotted lines) and two-body effective interactions derived 
from the N^LO NN potential (full lines). The Coulomb potential is included. 




FIG. 2: (Color online) The same as in Fig. 0for the point-nucleon radius. 
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FIG. 3: (Color online) Basis size dependence of the *He point-nucleon radius obtained using the two-body effective interactions 
derived from the N'^LO NN potential. Results for the harmonic-oscillator frequencies hfl = 24, 28, 32, 36, 40 MeV are presented. 



12 



'Li Minnesota 



10 11 12 13 14 15 16 

ha [MeV] 



17 



18 



- 
-A— 2 




19 



20 



21 



FIG. 4: (Color online) Harmonic-oscillator frequency dependence of the ®Li ground-state energy obtained in OTiQ-lAhQ basis 
spaces using two-body effective interactions derived from the Minnesota NN potential. Coulomb potential is not included. The 
dotted line corresponds to the Stochastic Variational Method result 
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FIG. 5: (Color online) Harmonic-oscillator frequency dependence of the ®Li ground-state energy obtained in 2hfl-14:hfl and 
8?ir2-16?if2 basis spaces using two-body effective interactions derived from the N"^LO NN potential (full lines) and the bare 
N^LO NN potential (dotted lines), respectively. The Coulomb potential is included. 
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FIG. 6: (Color online) Basis size dependence of the ®Li ground-state energy obtained in OTifl-lAfiQ and Ohfl-lGfiQ, basis spaces 
using two-body effective interactions derived from the N^LO NN potential and the bare N^LO NN potential, respectively. The 
Coulomb potential is included. The two-body effective interaction results correspond to the ground-state energy minimum for 
each model space. For details on the extrapolation see the text. 
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FIG. 7: (Color online) Basis size dependence of the ®Li ground-state energy obtained in Ohfl-lAhQ {16hQ, for MN) basis spaces 
using two-body effective interactions derived from the AV8', N^LO, CD-Bonn 2000 and the Minnesota NN potentials. The 
AV8' result obtained by the GFMC method [l^ and the Minnesota results obtained by the SVM ^ and the EIHH methods 
[sH are shown for a comparison. 
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FIG. 8: (Color online) Calculated positive-parity excitation spectra of ®Li obtained in 0hQ,-14fiQ, {12hQ for the 2^1 and ijo 
states) basis spaces using two-body effective interactions derived from the N'^LO NN potential are compared to experiment. 
The HO frequency of TiQ, = 8 MeV was used. The experimental values are from Ref. |39|. 
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FIG. 9: (Color online) The same as in Fig. IHlfor the HO frequency of hO. = 10 MeV. 
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FIG. 10: (Color online) The same as in Fig. El for the HO frequency of HQ. = 13 MeV. 
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FIG. 11: (Color online) Calculated positive-parity excitation spectra of ^Li obtained in the UTlQ {12hQ for the 2+1 and ijo 
states) basis space using two-body effective interactions derived from the N^LO NN potential are compared to experiment. 
The HO frequency dependence in the range of HQ = 8 — 13 MeV is presented. The experimental values are from Ref. )39| . 
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FIG. 13: (Color online) The dependence of the 1+0 and 3+0 state energy on the HO frequency for the 4hil, 6hQ and Sftil 
model spaces calculated using the N'^LO NN potential. In the inset, the 1+0 and 3+0 state energies at their respective HO 
frequency minima are plotted as a function of A'max- 
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FIG. 14: (Color online) Calculated positive-parity excitation spectra of ^"B obtained in ORfl-Sftn basis spaces using two-body 
effective interactions derived from the N^LO NN potential are compared to experiment. The HO frequency of h£l = 12 MeV 
where the 87if2 ground state is at the minimum was used. The experimental values are from Ref. [ssll . 
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Exp 


N^LO 


|£gs 1(1 + 0) 


31.995 


28.5(5) 


Qgs [e fm^] 


-0.082(2) 


-0.08(2) 




-hO.822 


-^0.839(5) 


£,(1+0) 


0.0 


0.0 


£;x(3+0) 


2.186 


2.91(3) 


Sx(0+1) 


3.563 


3.30(10) 


i:x(2+0) 


4.312 


4.10(15) 


Sx(2+1) 


5.366 


5.50(30) 


Sx(lJO) 


5.65 


5.40(30) 


B(E2;l+0 ^ 3+0) 


21.8(4.8) 


16.0(1.5) 


B(M1;0+1 ^ 1+0) 


15.42(32) 


15.01(10) 


B(E2;2+0 1^0) 


4.41(2.27) 


6.2(8) 


"He 


Exp 


N^LO 


|£gs 1(0+1) 


29.269 


26.2(5) 


"He^^Li 


Exp 


N^LO 


B(GT;0+1 1+0) 


4.728(15) 


5.22(10) 



TABLE II: Experimental and calculated energies, in MeV, quadrupole and magnetic moments, as well as E2, in fm**, and 
Ml, in fi%, transitions of "Li and "He as well as the B(GT) values for the "He ground state to "Li ground state transition. 
Results obtained using the N'^LO NN potential are presented. The errors are estimated from the HO frequency and basis size 
dependences as well as the bare interaction result extrapolations. The experimental values are from Ref. j3flll4nl|. 
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Exp 


N^LO 


CD-Bonn 


AV8' AV8' 


MN MN MN 


ISgeKl 


+0) 


31.995 


28.5(5) 


29.35(40) 


28.5(5) 28.19(5) 


36.50(4) 36.51 36.64(7) 


£;x(i+ 


0) 


0.0 


0.0 


0.0 


0.0 0.0 




Sx(3+ 


0) 


2.186 


2.91(3) 


2.86 


2.97 3.21(7) 




£;x(o+ 


1) 


3.563 


3.30(10) 


3.29 


3.41 3.94(8) 




£x(2+ 


0) 


4.312 


4.10(15) 


4.47 


4.28 4.10(6) 





TABLE III: Experimental and calculated energies, in MeV, of Li. Binding and excitation state energy results for the realistic 
N^LO, CD-Bonn 2000, AV8' NN potentials and the binding energy for the semi-realistic Minnesota (without Coulomb) NN 
potential are presented. The GFMC [T^, the SVM ^ and the EIHH ^ results are shown for comparison. See the text for 
further details. 



lOg 


Exp 


N^LO 


|£(3+0)| 


64.751 


54.65 


|£(1+0)| 


64.033 


56.28 


-Ex (3+0) 


0.0 


0.0 


-Ex(l^O) 


0.718 


-1.628 


£;x(o+i) 


1.740 


0.293 


Bx(lJO) 


2.154 


1.459 


Sx(2+0) 


3.587 


1.356 



TABLE IV: Experimental and calculated energies, in MeV, of ^'^B. The NCSM results obtained using the two-body effective 
interaction derived from the N^LO NN potential in the 8?!,n basis space and the HO frequency of hfl = 12 MeV. The experimental 
values are from Ref. [39| 



